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We report the specific heat (SH) measurements on single crystals of hole doped FeAs-based superconductor 
Bao (,Ko.4Fe2As2. It is found that the electronic SH coefficient y e (T) is not temperature dependent and increases 
almost linearly with the magnetic field in low temperature region. These point to a fully gapped superconducting 
state. Surprisingly the sharp SH anomaly AC/T\ Tc reaches a value of 98 mJ/molK 2 suggesting a very high 
normal state quasiparticle density of states (y„ « 63mJ/molK 2 ). A detailed analysis reveals that the y e (T) 
cannot be fitted with a single gap of s-wave symmetry due to the presence of a hump in the middle temperature 
region. However, our data indicate that the dominant part of the superconducting condensate is induced by an 
s-wave gap with the magnitude of about 6 meV. 

PACS numbers: 74.20.Rp, 74.25.Bt, 65.40.Ba, 74.70.Dd 



The discovery of high temperature superconductivity in 
the FeAs-based system has stimulated enormous interests in 
the field of condensed matter physics and material sciences 
[Q]]. The superconductivity has not only been discovered 
in the electron doped samples, but also in the hole-doped 
onesjj, Lj. The central issues concerning the superconduc- 
tivity mechanism are about the symmetry and the magnitude 
of the superconducting gap. The experimental results ob- 
tained so far are, however, highly controversial. The low 
temperature specific heat (SH) measurements in the F-doped 
LaFeAsO samples revealed a nonlinear magnetic field depen- 
dence of the SH coefficient y e , which was attributed to the 
presence of a nodal gap[4]. This was later supported bymany 

other measurements based on jusrih s n, nmrub m, 

magnetic penetration[ll] and point contact Andreev spec- 
trum (PCAS)[12]. On the other hand, the PCAS on the F- 
doped SmFeAsO indicated a feature of s-wave gap[13], some 
measurements ffl [l5l fl6l [l7tl also gave support to this con- 
clusion. It is important to note that most of the conclusions 
drawn for a nodal gap were obtained on the electron doped 
LnFeAsO samples (abbreviated as FeAs-\\\\, Ln stands for 
the rare earth elements) which are characterized by a low 
charge carrier density and thus low superfluid density! 18]. 
For the FeAs-\\\\ phase, it is very difficult to grow crys- 
tals with large sizes, therefore most of the measurements on 
the pairing symmetry so far were made on polycrystalline 
samples. This is much improved in the (Ba, S r)\- x K x Fe2A s 2 
(denoted as FeAs-122) system since sizable crystals can be 
achieved lll9[ 12011 . Preliminary data by angle resolved photo- 
emission spectroscopy (ARPES) on these crystals show two 
groups of superconductin g gaps ( Aj w 12 meV, A2 ~ 6 meV) 
all with s-wave symmetry l2l[ 22 . 2311 . It is known that the sur- 
face of this type of single crystals decay or reconstruct very 
quickly, this may give obstacles to get repeatable data when 
using the surface sensitive tools. Thus solid conclusions about 
the gap symmetry and magnitude from bulk measurements are 
highly desired. 

Specific heat (SH) is one of the powerful tools to mea- 
sure the quasiparticle density of states (DOS) at the Fermi 



level. By measuring the variation of the electronic SH ver- 
sus temperature and magnetic field, one can essentially deter- 
mine the feature of the gap symmetry. In this Letter, for the 
first time, we report the detailed low temperature SH data on 
BaudK()AFe2As 2 single crystals with T c = 36.5 K (90%p„). 
Our results elucidate a fully gaped feature of the supercon- 
ducting state. Meanwhile we show the evidence of a large 
DOS for the normal state of the FeAs-122 superconductors, 
which is in sharp contrast with that of the FeAs-l 1 1 1 phase. 

The superconducting single crystals with T c of about 36.5 
K were grown by the self-flux mefhod[20]. The sample for 
the SH measurement has the dimensions of 3 x 1.5 x 0.2 
mm 3 . The resistivity and the specific heat were measured with 
a Quantum Design instrument physical property measurement 
system (PPMS) with the temperature down to 1.8 K and the 
magnetic field up to 9 T. We employed the thermal relaxation 
technique to perform the specific heat measurements. To im- 
prove the resolution, we used a latest developed SH measuring 
puck from Quantum Design, which has negligible field depen- 
dence of the sensor of the thermometer on the chip as well as 
the thermal conductance of the thermal linking wires. 

In the main panel of Fig. 1 we show the raw data of SH 
coefficient y = C/T vs T at T and 9 T. Multiple compli- 
cated contributions to the SH data emerged in the low-T re- 
gion when a magnetic field was applied, so we only showed 
the data above 4.3 K under magnetic fields, and the data at 
zero field was shown down to about 1.8 K. Clear and sharp 
superconducting anomalies can be seen near T c from the raw 
data. The SH anomaly AC/T\t c at zero field was determined 
to be about 98 mJ/mol K 2 , indicated by the vertical short blue 
line in the inset (a). This is remarkably different from the case 
in the FeAs-l 1 1 1 phase, where no visible or o nly very small 
SH anomaly were observed in the raw datafl [24 . 251. Even 
in the same FeAs-122 system, the magnitude of the anomaly 
in our data is several times larger than that observed in the 
polycrystalline samples j26ll and the single crystals grown us- 
ing Sn as the flux ! 1 911 . This large value of AC/T\t c clearly 
suggests a rather high normal state quasiparticle DOS in this 
system, which will be further addressed later. A magnetic 
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FIG. 1: (color online) Raw data of SH coefficient y = C/T vs T are 
shown in the main frame. The dashed line shows the normal state 
SH obtained from fitting to eq.(4). The inset (a) shows an enlarged 
view of the data y = C/T near T c . The sharp SH anomaly AC/T\ Tc 
is indicated by the arrowed blue short line with a magnitude of about 
98 mj/mol K 2 . The inset (b) presents the resistive transition curves 
at magnetic fields ranging from to 9 T. 

field of 9 T shifts the SH anomaly down for only 1 .5 K and 
suppresses the anomaly. In the low-T region, a clear flatten- 
ing feature of C/T can be seen in the zero field data, which 
may imply the weak excitation of quasiparticles. The temper- 
ature dependence of the resistivity at different magnetic fields 
are shown in the inset (b) of Fig. 1. By applying a mag- 
netic field the middle transition point (50%p„) shifts to lower 
temperatures slowly with a slope -dp,oH C 2(T)/dT\T c ~ 4.1 
T / K. Using the Werthamer-Helfand-Hohenberg relation l27ll 
poH C 2(0) = -0.69dpoH C 2(T) / dT\rT c , we get the upper criti- 
cal field fi H c2 (0) * 100 T (H\\c). 

The raw data in the low-T region at different fields are plot- 
ted as C/T vs T 2 in Fig. 2(a). Detailed analysis reveals that 
weak Schottky anomaly still contributes to the whole SH in 
the low-T region. A slight curvature was also detected in the 
region from 4.3 K to 1 1 K on the plot of C/T vs. T 2 , which 
was attributed to the electron contribution of the supercon- 
ducting state and the quintic term of the phonon contribution. 
This brought in enormous difficulties when treating the data 
because it gave too many fitting parameters. Therefore we 
first analyze the data below 6 K at zero field, where the two 
terms mentioned above remain negligible. Consequently the 
zero field data below 6 K can be represented by the following 
equation: 

C(T,H = 0) = y T +/3T 3 + C Sch (T,H= 0), (1) 

where the three terms represent the contributions of the resid- 
ual electronic SH, the phonon and the magnetic impurity 
(the so-called Schottky anomaly), respectively. The two- 
level Schottky anomaly is given by nx 2 e x /(l + e x ) 2 (x = 
gHBH e ff/k B T), where g is the Lande factor, fi B is the Bohr 

magneton, H e //= JH 2 + H 2 is the effective magnetic field 
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FIG. 2: (color online) Temperature and magnetic field dependence of 
specific heat in C/T vs T 2 plot in the low-T range, (a) Raw data be- 
fore removing the Schottky anomaly. The dashed lines represent the 
theoretical fit (see text) containing all terms in Eq. (2). (b) Replot of 
the data after the Schottky anomaly was subtracted. The dot-dashed 
line represents a extension of the zero field data to T = K giving a 
residual value y = 7.7 ml/mol K 2 (see text). 

which evolves into H e ff - Ho, the crystal field at zero external 
field, and n is the concentration of paramagnetic centers. The 
obtained fitting parameters y ( > as 7.7 mJ/mol K 2 and /3 ~ 0.473 
mJ/mol K 4 are very close to the values obtained by simply 
drawing a linear line below 6 K as shown by the dot-dashed 
line in Fig. 2(a). The values of jq and /3 are then fixed when 
fitting the zero field data up to 11 K, where all the terms must 
be taken into account: 

C(T, H = 0) = y T + IJ3T 3 +j 1 T 5 ]+C es + C S ch(T, H = 0), (2) 

where rj is the quintic term coefficient of the phonon SH and 
C es — D x e~ A( - ^ k " T /T L5 is the superconducting electron con- 
tribution, with A(0) the superconducting gap at K. By fitting 
the data at zero field using Eq. (2), we obtained 77 « 0.00034 
mJ/mol K 6 and A(0) ~ 5.99 ± 0.3 meV. As for the data under 
finite fields, a magnetic field induced term y(H) arises and the 
total SH can be written as 

C(T,H) = [y +y(H)]T + [/3T 3 +j ] T 5 ]+C es + C Sch (T,H). (3) 

It is quite rational to fix yo, /3, and 77 as the values obtained 
from analyzing the data at zero field. And the supercon- 
ducting gap under magnetic fields was restricted using the 
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FIG. 3: (color online) Field dependence of the field-induced term 
y{H) at temperatures ranging from 4.5 K to 11.1 K, and that at T = 
K obtained from fitting (see text). The dashed lines for different 
temperatures are guides for the eyes. The blue solid line is a linear 
fit to the zero temperature data, and the red dotted line is a fit to the 
d-wave prediction y(H) = A \Th. 



FIG. 4: (color online) (a) Scaling of the data according to the s-wave 
scenario (symbols) C cai - S vs. T/ \fH, the dashed line represents the 
theoretical expression, (b) Scaling of the data (symbols) based on 
d-wave prediction C ca i-a vs. T/ Vff. No good scaling can be found 
for the d-wave case. 



TABLE I: Fitting parameters (A is calculated through A(H) = 
A(0) VI - H/H C 2, Ho is fixed with the value from the fitting to zero- 
field data). 

/4)ff(T) O(mJK - 5 /mol) A(meV) n(mJ/molK) fi H (T) g 



0.0 


2.11 x 10 6 


5.99 


20.95 


1.70 


2.75 


1.0 


2.07 x 10 6 


5.96 


31.08 


1.70 


3.18 


2.0 


1.93 x 10 6 


5.93 


35.76 


1.70 


3.26 


3.0 


1.82 x 10 6 


5.90 


36.22 


1.70 


3.38 


5.0 


1.63 x 10 6 


5.84 


36.11 


1.70 


3.24 


7.0 


1.57 x 10 6 


5.78 


36.71 


1.70 


3.11 


9.0 


1.43 x 10 6 


5.72 


35.98 


1.70 


3.22 



relation|28J A(H) = A(0) \A - H/H c2 assuming a field in- 
duced pair breaking effect. In this way the number of the 
fitting parameters were reduced remarkably and creditable re- 
sults can be obtained. 

Fig. 2(b) shows the data after removing the Schottky 
anomaly from the total SH. The obtained field induced term 
y(H) are shown in Fig. 3 (will be discussed later) and the fit- 
ting parameters related to the terms C es and Csch are shown 
in Table I. The obtained residual term yo ~ 7.7 mJ/mol K 2 
accounts for about 1 1 % of the total electron contribution (will 
be discussed later), indicating a superconducting volume frac- 
tion of about 89% in our sample. Using the obtained value 
of p and the relation & D = (l2n 4 k B N A Z/5p) 1/3 , where N A = 
6.02 xlO 23 moF 1 is the Avogadro constant, Z = 5 is the num- 
ber of atoms in one unit cell, we get the Debye temperature 
®o ~ 274 K. This value is comparable to that found in the 
LaFeAsO0.9F0.1_i system^. 

The field-induced change of the electron SH coefficient 
y(H) was investigated carefully. This term at zero tempera- 
ture obtained from fitting, along with the data at finite temper- 



atures extracted from Fig. 2(b), is plotted in Fig. 3. It can be 
seen clearly that y(H) increases almost linearly with the mag- 
netic field in the temperature region up to 1 1 K. A linear fit 
with the slope of about 0.633 mJ/molK 2 T to the zero tem- 
perature data is revealed by the blue solid line in this figure. 
This linear behavior is actually anticipated by the theoretical 
prediction for superconductors with a full gapj29j], in which 
y{H) is mainly contributed by the localized quasiparticle DOS 
within vortex cores. This is in sharp contrast with the results 
in cuprates[30, 31] and the LaFeAsO0.9F0.1-5 system[4] where 
a y(H) oc y[H relation was observed and attributed to the 
Doppler shift of the nodal quasiparticle spectrum. The curve 
plotted using the relation y(H) - A yflJ is also presented in 
Fig. 3 by the red dotted line for comparison. It is obvious that 
this curve fails to fit our data. 

In order to further confirm the gap symmetry, we ana- 
lyzed the SH data in finite temperature region in the mixed 
state. In s-wave superconductors, the inner-core states domi- 
nate the quasiparticle excitations, and consequently a simple 
scaling law C core /T 3 * (y„/7f c2 (0)) x (T/ \/#T 2 is expected 
at low-T. While for a gap with line nodes, the excitation spec- 
trum is dominated by the extended quasiparticles outside the 
vortex cores. And the so-called Simon-Lee scaling law[32] 
C vo i/(T yfH) = f(T/ V77) may be obeyed. A simple anal- 
ysis similar to that has been done in our previous work lO 
and shows that for the superconductor with an s-wave sym- 
metry, C C ais = [(C(H) - Csch(H)) - (C(H = 0) - C Sc h(H = 
0))]/r 3 ss C core /T 3 . In other words, the defined term C ca \- S 
should scale with (77 VTfT 2 with the prefactor y w /77 C 2(0). 



33] 



Similarly for the d-wave symmetry we have known[31 
that Cca-d = [(C(H) - C Sch {H)) - (C(H = 0) - C Sch (H = 
0))]/rV/7 = C vo i/TVH - aT/y/H, where a is the elec- 
tron SH coefficient at zero field for a d-wave superconductor, 
should also scale with T/ y/H. 
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FIG. 5: (color online) Temperature dependence of the electronic SH 
contribution (with the normal state part subtracted) is shown in the 
main frame. A sharp SH anomaly can be seen here. A hump is 
clearly seen in the middle temperature region. The red dashed line is 
a theoretical curve based on the BCS expression with an s-wave gap 
of 6 meV. The inset shows the entropy of the superconducting state 
(red circle symbols) and the normal state (dark square symbols). 



The scaling result with the s-wave condition is presented in 
Fig. 4(a). One can see that all the data at different magnetic 
fields can be scaled roughly to one straight line, which reflects 
the theoretical curve C cal - S = 0.633 x (T/ V#r 2 . Naturally, 
this prefactor y„/H C 2(0) = 0.633 mJ/mol K 2 T is consistent 
with the magnitude of the slope of the blue line in Fig. 3. 
Using the value of H C 2(0) ~ 100 T, we can estimate the nor- 
mal state electron SH coefficient y„ of about 63.3 mJ/mol K 2 . 
Fig.4(b) shows the scaling by following the d-wave scheme. 
It is clear that the s-wave scaling is much better than that of 
the d-wave case. So it seems that the superconducting gap in 
this FeAs-122 phase has an s-wave symmetry and the field- 
induced quasiparticle DOS are mainly contributed by the vor- 
tex cores. 

Using the value of y„ * 63.3 mJ/mol K 2 , we get the ra- 
tio AC e /y n T\T c ~ 1-55 being very close to the weak-coupling 
BCS value 1.43. Considering the electron re-normalization 
effect, the electron SH coefficient of a metal can be written as 
y„ = 2n 2 N{E F )k 2 B {\ + A)/3, where N(E F ) is the DOS at the 
Fermi surface and A reflects the coupling strength. The fact 
that the electron-phonon coupling strength is weak in present 
system indicates that the large value of y„ is not due to the 
enhanced effective mass but originates from the high normal 
state quasiparticle DOS. Comparing with the y„ obtained for 
the F-doped LaFeAsO system] 3l (about 5-6 mJ/mol - FeK 1 ), 
the N(Ef) in hole doped FeAs-122 may be 3-5 times higher 
than that in the electron doped FeAs-1 1 11. This may give an 
important clue in the calculation of the N(Ef) induced by the 
doping effect. 

In the FeAs-122 superconductors, it is challenging to mea- 
sure the normal state SH below T c due to the very high H C 2- 
In order to have a comprehensive understanding to the normal 
state electronic SH, we have attempted to fit the normal state 



SH above T c using a polynomial function: 

C„ = (y +y„) +p 3 T 3 +/3 5 T 5 +y 6 7 r 7 +/3 9 T 9 +0 u T n , (4) 

where we took the values obtained already yo = 7.7 mJ/mol 
K 2 , = 0.473 mJ/mol K 4 , y„ = 63.3 mJ/mol K 2 . Other fit- 
ting parameters, {I5, and /3\ \ , were left free in the fitting 
process, yielding the values of 3.72x 10~ 4 mJ/mol K 6 , -5.32x 
10~ 7 mJ/mol K 8 , 2.13X 10" 10 mJ/mol K 10 , and -2.90x 10~ 14 
mJ/mol K 12 , respectively. It's worth to note that the value of 
/?5 is very close to the value of rj obtained before. The fitting 
result of the normal state SH is displayed by the blue dashed 
line in the main frame of Fig. 1 . The data after subtracting the 
normal state SH is presented in the main frame of Fig. 5. It 
was found that the entropy-conserving law was satisfied nat- 
urally confirming the validity of our fitting, as shown in the 
inset of Fig. 5. A clear flattening of y e - y„ in the temper- 
ature region up to 7 K is observed indicating a fully gapped 
superconducting state. Moreover, a hump is clearly seen in 
the middle temperature region. We attempted to fit the data 
using the BCS formula: 



4jV(0) 
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where f = y/s 2 + A 2 (T, 0) and A(T, 0) = A Q (T) for the s-wave 
symmetry. The red dashed line in the main frame presented 
the fitting result. One can see that the fitting curve with a gap 
value of about 6 me V matched our data below 1 3 K perfectly, 
but failed to describe the hump feature in the middle temper- 
ature region. This hump may be attributed to the multi-gap 
effect which seems to appear in the FeAs based supercon- 
ductors. But at this moment we can't exclude the possibility 
that the hump is induced by the limited uncertainty in getting 
the normal state phonon contribution. Nevertheless, the fine 
fitting in wide temperature region strongly suggests that the 
dominant part of the superconducting condensate is induced 
by an s-wave gap with the magnitude of about 6 meV. Our 
results here seem to be consistent with the ARPES data, both 
in symmetry and the small gap j2lll22ll23ll . But we have not 
found a large gap of 12 meV. This discrepancy may be induced 
by the different ways in determining the gap. Future works are 
certainly required to reconcile all these distinct results. 

In summary, the low temperature specific heat measure- 
ments reveal that the Bao.6^0.4^2^2 superconductor has an 
s-wave pairing symmetry with the gap amplitude of about 
6 meV. The rather high SH anomaly AC/T\t c and the large 
value of y„ suggest a high normal state quasiparticle DOS in 
this system. This makes it very different from the FeAs-1 1 1 1 
phase. A multigap feature seems possible but the other gap 
should be responsible for only a small fraction of the super- 
fluid density. 
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